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The dynamic magnetization-reversal phenomena in the Ising model under a finite-duration exter- 
nal magnetic field competing with the existing order for T < T° has been discussed. The nature of 
the phase boundary has been estimated from the mean-field equation of motion. The susceptibility 
and relaxation time diverge at the MF phase boundary. A Monte Carlo study also shows diver- 
gence of relaxation time around the phase boundary. Fluctuation of order parameter also diverge 
near the phase boundary. The behavior of the fourth order cumulant shows two distinct behavior: 
for low temperature and pulse duration region of the phase boundary the value of the cumulant 
O^l , at the crossing point for different system sizes is much less than that corersponding to the static 

transition in the same dimension which indicate a new universality class for the dynamic transition. 
Also, for higher temperature and pulse duration, the transition seem to fall in a mean-field like 
weak-singularity universality class. 
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Thin films of ferromagnetic materials are of considerable interest as potential materials for the use as ultra-high 
' density recording medium. Recently, sophisticated experimental techniques such as the magnetic force microscopy 
have enabled physicists to look at the magnetization state and switching mechanism of ferromagnetic particles in the 
nanometre scale Q . Experiments 0, Q on i ron sesquilayers grown on tungsten are interesting from the theoretical 
point of view as the results seem to correspond precisely to an ideal ferromagnetic Ising system Q . 

The response of a pure magnetic system to time-dependent external magnetic fields has been of current interest in 
| statistical physics 0,1HIS0 These studies, having close applications in recording and switching industry, have also 
^ ■ got considerable scientific and technical interest. These spin systems, driven by time-dependent external magnetic 
fields, have basically got a competition between two time scales: the time-period of the driving field and the relaxation 
time of the driven system. This gives rise to interesting non-equilibrium phenomena. Tome and Oliveira first made a 
mean-field study of kinetic Ising systems under oscillating field. The existence of the dynamic phase transition for 
such a system and its nature have been thoroughly studied using extensive Monte-Carlo simulations [a, Q • Recently, 
• Junier and Kurchan [ljj have shown that the dynamics under oscillating field can induce asymmetry in the order 
parameter; namely an isotropic Heisenberg model (classical) under an alternating field switches to XY-like or Ising-like 
"p-; behavior, as the external field frequency and amplitude is tuned. 

Later, investigations were extended to the dynamic response of (ferromagnetic) pure Ising systems under magnetic 
fields of finite-time duration 0, 0. All the studies with pulsed field were made below T°, the static critical 
temperature (without any field), where the system gets ordered. A 'positive' pulse is one which is applied along the 
direction of prevalent order, while the 'negative' one is applied opposite to that. The results for the positive pulse 
case did not involve any new thermodynamic scale [TT| . In the negative pulse case, however, interesting features 
q [ were observed ^2 HS> ITU fl6| : the negative field pulse competes with the existing order, and the system 
• • ■ makes a transition from one ordered state characterised by an equilibrium magnetisation +mo (say) to the other 
equivalent ordered state with equilibrium magnetisation —mo, depending on the temperature T, field strength h p and 
its duration At. This transition is well studied in the limit At — ► oo for any non-zero value of h p at any T < T c °. This 
transition, for the general cases of finite At, is called here 'magnetisation-reversal' transition. As in the oscillating 
field case, in the study of magnetizatin-reversal transition, the finite duration (At) of the external field introduces 
a competition between two time scales: the relaxation time of the system and At. Some aspects of this transition 
has recently been studied extensively [f| 0, El 0] • Recently, studies on switching behavior of some thin films of 
'perpendicular media' (where spins are Heisenberg like rather than Ising like) have identified the same transition and 
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FIG. 1. Schematic diagram showing the free energy landscape against the magnetization at different times. For t < to the landscape 
is symmetric with respect to m = and the system is in +mo state as shown by • . For to < f < to + At, the landscape acquires a 
global minimum at m = —mo due to the presence of the external field and the system starts moving towards it. Beyond t = to + At, the 
landscape again becomes symmetric and the system eventually stabilizes in either of the equivalent minima as shown by o. 



II. MODEL AND THE TRANSITION 



The model studied here is the Ising model with nearest-neighbour interaction under a time-dependent external 
magnetic field. This is described by the Hamiltonian: 

H=-j2 J iJ s i s 3- h (t)Yl s i> w 

{y} « 

where Jy is the cooperative interaction between the spins at site i and j respectively, and each nearest-neighbour pair 
is denoted by {...}. We consider a square lattice. The static critical temperature is Tj? = 2/ln(l + v2) — 2.269... (in 
units of J I Kb)- At T < T c °, an external field pulse is applied, after the system is brought to equilibrium characterised 
by an equilibrium magnetisation mo(T). The spatially uniform field has a time-dependence as follows: 

h{t) = \-hv t <t<t + At 
w [0 otherwise. ^ ' 




FIG. 2. Typical time variation of the response magnetisation m(t) for two different field pulses h(t) with same At for an Ising system 
at a fixed temperature T. The magnetisation-reversal here occurs due to increased pulse strength, keeping their width At same. The 
transition can also be brought about by increasing At, keeping h p fixed. The inset indicates the typical phase boundaries (where the field 
withdrawal-time magnetisation m w = 0) for two different temperatures (sequential updating; note that for random updating the phase 
boundaries shift upwards). 
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Typical time-dependent (response) magnetisation m(t) (— < Si >, where < ... > denotes the thermodynamic 
'ensemble' average) of the system under different magnetic field h(t) are indicated in the Fig. 2. The time to at which 
the pulse is applied is chosen such that the system reaches its equilibrium at T (< T c ° ). As soon as the field is applied, 
the magnetisation m(t) starts decreasing, continues until time t + At when the field is withdrawn. The system relaxes 
eventually to one of the two equlibrium states (with magnetisation —too or +toq). At a particular temperature T, 
for appropriate combinations of h p and At, a magnetisation-reversal transition occurs, when the magnetisation of the 
system switches from one state of equilibrium magnetisation too to the other with magnetisation —too. This reversal 
phenomena at T < Tj? is simple and well studied for At — > oo for any non-zero h p . We study here the dynamics 
for finite At values, ft appears that generally h p — > oo as At — > and h p — > as At — > oo for any such dynamic 
magnetisation-reversal transition phase boundary at any temperature T (< T®). In fact, a simple application of the 
domain nucleation theory gives h p \nAt = constant along the phase boundary, where the constant changes by a factor 
l/(d-\- 1), where d denotes the lattice dimension, as the boundary changes from single to multi-domain region |ll| . 

III. MEAN-FIELD STUDY 

Consider a system of N Ising spins in contact with a heat bath evolving under Glauber single spin flip dynamics. 
The master equation can be written as |l8j 



-P (Si, ■■■,S N ;t) = -J2 W j (S j )P (Si, ...,S N ;t) + J2 W 3 (-S 3 )P (Si, .... S, S N ; t) , 



where P (Si, • • • , Sn; t) is the probability to find the spins in the configuration (Si, • • • , Sn) at time t and Wj(Sj) is 
the probability of flipping of the jth spin. Satisfying the condition of detailed balance one can write the transition 
probability as 



Wj(Sj) = & [l - Sj tanh ( S«W>+*» 
where A is a temperature dependent constant. Defining the spin expectation value as 



m% = (Si) = Y,s s * p ( s i> ■ ■ ■ ' t) , 
where the summation is carried over all possible spin configurations, one can write 



drrii 

~df = ~ m > 



ytanh 



12 j JijSj + hi 
f 



(3) 



Using the mean field approximation linearisation (in to), (PJ can be written after a Fourier transform, 



dm q (t) 
dt 



= A" 



(K(q)-l)m q (t) 



h g (t) 



T 



(4) 



where K(q) — J(q)/T. When we are concerned only with the homogeneous magnetization, we consider the q = 
mode of the equation and writing to 9= q = m an d h q= o = h, we get 



dm 
~dt ~ 



(5) 



In the mean field approximation K(0) = T^ IF /T with T C MF = J(0) and for small q, K(q) ~ K(0) (l - q 2 ). Differen- 
tiating Q with respect to the external field, we get the rate equation for the dynamic susceptibility x 9 (t) as 



d Xg{t) = A _i 
dt 



(K(q)-l) Xq (t) + - 



(6) 
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We divide the entire time zone in three different regimes : (I) < t < to, where h(t) = (II) to < t < to + At, 
where h(t) — — h p and (III) t + At < t < oo, where h(t) — again. We note that 10 can be readily solved separately 
for the three regions as the boundary conditions are exactly known. In region I, dra/dt — and the solution of the 
linearized 10 becomes trivial. We, therefore, use the solution of (QJ in region I (mo = tanh (moT^ F /T)) as the initial 
value of m for region II. Integrating in region II, we then get 

m{t) = % + ( m ° " at) exp [6AT (t " to)] ' (7) 

where b = 1/AT and AT = T C MF — T. We also note that in order to justify the validity of the linearization of J3J one 
must keep the factor inside the exponential of (7J small. This restricts the linear theory to be valid close to T^ IF and 
for small At. Writing m w = m(t + At), we get from Q 



h p ( hp \ bATAt , , 

~~ AT V _ AT I ^ ' 

It is to be noted here that in absence of fluctuations, the sign of m w (h p , At) solely decides which of the final two 
equilibrium states will be chosen by the system after the withdrawal of the pulse. At t = t + At, if m w > 0, the 
system gets back to +mo state and if m w < 0, magnetization- reversal transition occurs and the system eventually 
chooses the —mo state (see figure 2). Thus setting m w — 0, we obtain the threshold value of the pulse strength at the 
mean field phase boundary for this dynamic phase transition. At any T, combinations of h p and At below the phase 
boundary cannot induce the magnetization-reversal transition, while those above it can induce the transition. From 
(JHl therefore we can write the equation of the mean field phase boundary for the magnetization-reversal transition as 




1 10 100 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 

At m w 



FIG. 3. (a) MF phase boundaries for three different temperatures. The solid line is o btain ed from numerical solution of in MF case and 
the dotted lines give the corresponding analytical estimates in the linear limit (from Il4l ). (b) Logarithmic divergence of t^ f across the 
phase boundary for T/T c = 0.9. The data points (circles) for t^ f are obtained from the s olut ion of (with 5 = 0, \m(t) — m w \ < 10~ 4 ) 
and the solid line corresponds to 1111 is the solution of the linearized MF equation (from |14| ). 

Figure 3a shows phase boundaries at different T obtained from © and compares those to the phase boundaries 
obtained from the numerical solution of the full dynamical equation (0J . The phase boundaries obtained under linear 
approximation match quite well with those obtained numerically for small values of At and at temperatures close 
to T^ F , which is the domain of validity of the linearized theory as discussed before. In region III, we again have 
h(t) = and solution of (J3J) leads to 



m(t) = m w exp [6AT {t - (t + At)}} 



(10) 
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We define the relaxation time Tjf F , measured from t = to + At, as the time required to reach the final equilibrium 
state characterized by magnetization ±m in region III (see figure 1). From above expression therefore we can write 



MF 



1 , I mo 
■ in 



b AT \ \m„ 
( T 



y T MF _ T 



In 



(11) 



A point to note is that m(t) in (|10() grows exponentially with t and therefore in order to confine ourselves to the 



linear regime of m(t), mo must be small (T close to T c ) and t < Tr . The factor (T ( 



-•MF 



critical slowing down for the static transition at T = T ( 



MF 



However, even for T <C T f 



MF 



r MF 

R 

r MF 



gives the usual 
diverges at the 



magnetization- reversal phase boundary where m w vanishes. Figure 3b shows the divergence of t^ * against m w as 
obtained from the numerical solution of the full mean field equation of motion 10} and compares it with that obtained 
from JTTJ- 

Solution of Xq(t) is more difficult as all the boundary conditions are not directly known. However, XgW can b e 
expressed in terms of m(t) and the solution of the resulting equation will then have the t dependence coming through 
m(t), which we have solved already. Dividing |HJ by (J5J) we get 



dm 



which can be rewritten in the linear limit as 



T(K(q)-l) 



*(*)( 



K(0)-1 
K(q)-1 



h(t) 
K(q)-1 



(12) 



dXq 



dm 



X q +T q m + h(t)/At' 

where V = 1/T (K(q) - 1) and a q = (K(q) - 1) / (if (0) - 1) ~ 1 - q 2 / AT for small q. 
In region II, solution of (fH5|> can be written as 



(13) 



Xq{t) 



(x s q + r) 



m(t) - h p / AT 



tuq — hp/ AT 



(14) 



where x q is the equilibrium value of susceptibility in region I. Solving (|14fl in region III with the initial boundary 
condition m (to + At) = m w , we get 



Xq (t) = -r + (x q (t + At) + r) 

' m(t) 



-r + ( Xq + r) 



m(t) 



"a 



^bATAta„ 



(15) 



where use has been made of (|14|) and JSJl. The dominating q dependence in Xq{t) IS coming from {l/m w ) aq when 
m w — > as one approaches the phase boundary. The singular part of the dynamic susceptibility can then be written 
as 



2 (<-MF\ 



Xq(t) = (^ + r)exp [-q 2 (£ 
where for small values of m w the correlation length £ MF is given by 0, 0| 



(16) 



AT \ m„ 



(17) 
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Thus the length scale also diverges at the magnetization-reversal phase boundary and this can be demonstrated even 
using the linearized mean field equation of motion. Equations IjllH and (|17|) can now be used to establish the following 
relation between the diverging time and length scales : 



.MF 



T (fMF\ 2 



(18) 



which leads to a dynamical critical exponent z — 2. One can have a numerical estimate of £ MF by solving (|16|1 for 
different values of q. Figure 4(b) shows plots of x g (i) against m w for different values of q. It may be noted that 
these divergences in tJj and £ MF are shown to occur for any T < T^ IF , and these dynamic relaxation time and 
correlation length defined for the magnetization-reversal transition exist only for T < T C MF . 

It may further be noted from (fT^fl that Xq(t) ~~ > as £ MF ~ * 00 > thereby producing a minimum of Xq at the phase 
boundary. The absence of any divergence in the susceptibility is due to the fact that at t = to + At, there remains 
no contribution of m w in x g (t) as ls evident from (| 1 5 p . However, numerical solution of (|12[1 for q — mode shows 
a clear singularity in the homogeneous susceptibility xo at the magnetization-reversal phase boundary (m w — 0), as 
depicted in figure 4a. 
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FIG. 4. (a) Divergence of Xq=0 across the phase boundary obtained from the numerical solution of 1121 . (b) Plot of \ q against m m 
for different values of q. The inset shows the linear variation of (£ MF ) 2 against [ln|muj|] — 1 . The data points for £ MF i n th e inset are 
obtained from the slope of the best fitted straight lines through a plot of ln\q against q 2 for different values of m w (from jlJI t. 



IV. MONTE-CARLO STUDY AND THE RESULTS 



Here the Monte-Carlo study has been carried out in 2D and 3D (square and cubic lattices) with periodic boundary 
conditions. Spins are updated following Glauber dynamics. The updating rule employed here are both sequential 
as well as random. In sequential updating rule one Monte-Carlo step consist of a complete scan of the lattice in a 
sequential manner; while in random updating a Monte-Carlo step is defined by N (— L d ) random updates on the 
lattice, where ./V is the total number of spins in a lattice of linear size L in d dimensions. Studies have been carried 
out at temperatures below the static critical temperature (e.g, T® ~ 2.27 for 2D). The system is allowed to evolve 
from an initial state of perfect order to its equilibrium state at temperature T. The time to is chosen to be much 
larger than the static relaxation time at that T, so that the system reaches an equilibrium state with magnetization 
-|-mo (T) before the external magnetic field is applied at time t = to- The field pulse of strength —h p is applied for 
duration At (measured in Monte Carlo steps or MCS). The magnetisation starts decreasing from its equilibrium value 
mo- The average value of the magnetisation m w at the time of withdrawal of pulse is noted. The phase boundary of 
this dynamic transition is defined by appropriate combination of h p and At that produces the magnetisation reversal 
by making m(to + At) = m w = from a value m(to) — too, i.e, m w changes sign across the phase boundary. The 
phase boundary changes with T. The behavior of different thermodynamic quantities are studied across the phase 
boundary. These quantities are averaged over different initial configurations of the system. The fluctuations over the 
average value are also noted. 

We note the marked qualitative difference between the MF and MC phase boundaries (see Fig. 3 and Fig. 5). In 
the MF case, even for At — > oo, due to absence of fluctuations, h p must be greater than the non-zero coercive field in 
order that the transition takes place and the phase boundary becomes fiat for large At. In MC study, due to presence 
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of fluctuations, even a pulse of infinitesimal strength, if applied for very long time, can induce the transition. This is 
evident from the asymptotic nature of the phase boundaries for large At. 
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FIG. 5. Phase boundaries from Monte Carlo studies on square lattice with L 
random updating (filled symbol). 
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To understand the nature of the MC phase diagram of the magnetization-reversal transition we consult the classical 
theory of nucleation (CNT). A typical of a ferromagnet, below its static critical temperature T°, consists of droplets 
or domains of spins oriented in the same direction, in a sea of oppositely oriented spins. According to CNT, the 
equilibrium number of droplets consisting of s spins is given by n s — A exp (— e s /T), where e s is the free energy of 
formation of a droplet containing s spins and N being the normalization constant. In presence of a negative external 



magnetic field h, the free energy can be written as e s 



-2hs 



as 



(d-l)/d 



, where the droplet is assumed to be 



spherical and c(T) is the temperature dependent surface tension. Droplets of size greater than a critical value s c are 
favoured to grow, where s c = [er(d — 1)/(2g? \h\)] d is obtained by maximizing e s . The number of supercritical droplets 



is therefore given by n Sc = N exp 



-k d a d \h\ 1 - d /T 



where A^ is a constant depending on dimension only. In the SD 

regime, where a single supercritical droplet grows to eat up the whole system, the nucleation time goes inversely as 
the nucleation rate /. According to the Becker-Doring theory, / is proportional to n Sa and therefore one can write 



t^ d (x I 1 oc exp 



T\h 



d-1 



However, in the MD regime the nucleation mechanism is different and in this regime many supercritical droplets grow 
simultaneously and eventually coalesce to create a system spanning droplet. The radius s^ d of a supercritical droplet 
grows linearly with time t and thus s c oc t d . For a steady rate of nucleation, the rate of change of magnetization is 
It d . For a finite change Am of the magnetization during the nucleation time tJ^ d , one can write 



It d dt = I (t 



MD\ d + 1 
N ) 



Therefore, in the MD regime one can write 19] 2 



l/(d+l) 



ex exp 



T(d+l)\h\ c 



(19) 



During to < t < to + At, when the external field is 'on', the only relevant time scale in the system is the nucleation 
time. The magnetization reversal phase boundary gives the threshold value hp of the pulse strength which, within 
time At, brings the system from an equilibrium state with magnetization +mo to a non-equilibrium state with 
magnetization m w — 0_, so that eventually the system evolves to the equilibrium state with magnetization —mo. 
The field driven nucleation mechanism takes place for to < t < to + At and therefore equating the above nucleation 
times with At, one gets for the magnetization-reversal phase boundary 



In (At) 



Cl 

C2 



l-d 



c[h; 



7(d- 



in the SD regime 
1), in the MD regime 



(20) 
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where C = Ad& d /T and c±, ci are constants. Therefore a plot of ln(At) against [h°] would show two different 
slopes corresponding to the two regimes ^J. Figure 6 shows these plots and it indeed have two distinct slopes for 
both d = 2 (figure 6(a)) and d — 3 (figure 6(c)) at sufficiently high temperatures, where both the regimes are present. 
The ratio R of the slopes corresponding to the two regimes has got values close to 3 for d = 2 and close to 4 for d = 3, 
as suggested by (|2(JH . The value of hjygp is obtained from the point of intersection of the straight lines fitted to the 
two regimes. At lower temperatures, however, the MD region is absent and the phase diagram here is marked by a 
single slope as shown in figures 6(b) and 6(d). 
The relaxation behaviour follows : 



t r ~ n(T,L)e 



(21) 



where k(T, L) is a constant depending on temperature and system size and n(T) is a constant depending on tem- 
perature only, ft may be noted from l|21|l that r — > k(T, L) as m w — > 0. Therefore the true divergence at the phase 
boundary (where m w = 0) of the relaxation time depends on the nature of k(T, L). k(T,L) grows sharply with 
system size. The relaxation time tr therefore diverges in the thermodynamic limit (L — * oo) through the constant 
k. It may be noted that this divergence of tr at the dynamic magnetization-reversal phase boundary occurs even at 
temperatures far below the static critical temperature T c °. 
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FIG. 6. Plots of In At against {h p ) 1 ~ d along the MC phase boundary, (a) T/T c = 0.31 and (b) T/T c = 0.09 for square lattice and (c) 
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T/Tc = 0.67 and (d) T/T c = 0.11 for simple cubic lattice. The slope ratio _R ~ 3.27 in (a) and ~ 3.97 in (c) (from Qj) 
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FIG. 7. Finite size scaling fits : (a) for d = 2 at T/T c = 0.88 and (b) for d = 3 at T/T c = 0.67 (from [Ti l. 

In the region where the transition is continuous, the scaling arguements are expected to hold. We assume power 
law behaviour in this regime both for m w 



\h p -h c p (At,T)\ 



and for the correlation length 



(22) 



^ \h p - h c (At,T)\ 



(23) 
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For a finite size system, hp is a function of the system size L. Assuming that at the phase boundary £ can at the 
most reach a value equal to L, one can write the finite size scaling form of m w as: 



m w ~ L-P/»f [(h p - h; (At, T, L)) L 1 ^] , (24) 

where f(x) ~ x^l v as x — ► oo. A plot of m w /L~Pl u against — /i£ (At, T, L)) L 1 /" shows a nice collapse of the data 
for d = 2 and for d = 3 as shown in figure 7. The values of the critical exponents obtained from the data collapse 
were (3 = 0.85 ± 0.05 and v = 1.5 ± 0.5 in d = 3 and /3 = 1.00 ± 0.05 and v = 2.0 ± 0.5 in d = 2, where (At, T) 
was obtained with an accuracy O (lO -3 ). Attempts to fit similar data to the above finite size scaling form obtained 
in the SD regime failed. 

The behavior of the reduced fourth order cumulant U |2jJ has also been studied near the magnetisation reversal 
transition boundary. This is defined as 

u = 1 _Jj<)_ ! (25) 
3 (ml) 

where (rn^) is the ensemble average of m^,. (jn^) is similarly defined. The cumulant U here behaves somewhat 
differently compared to that in static and other transitions: Deep inside the ordered phase m w ~ 1 and U — > 2/3. For 
other (say static) transitions the order parameter (m w ) goes to zero with a Gaussian fluctuation above the transition 
point, giving U — > there. Here, however, due to the presence of the pulsed field, \m w \ is non-zero on both sides of 
the magnetisation-reversal transition. Hence U drops to zero at a point near the transition and grows again after it. 





FIG. 8. Behavior of U near the transition, driven by (a) T at a fixed value of h p (=1.9) and At (=5) with sequential updating, (b) h p 
at a fixed value of T (=0.5) and At (=5) with sequential updating, and (c) h p at a fixed value of T (=0.5) and At (=5) with random 
updating, for different L, averaged over 1000 to 20000 initial configurations. The insets show the typical behavior of the magnetisation 
m w at the time of withdrawal of the field pulse by varying (a) T at a fixed h p and At, for L = 100 and 800, (b) h p at a fixed T and At, 
for L = 100 and 400, (c) h p at a fixed T and At, for L = 50 and 200; m w = at the effective transition point, (d) Finite size scaling 
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study in this parameter range: the effective T c or values (see the insets), where m w = 0, are plotted against L 1 / y with v 1 = 0.7. 
The values of the cumulant crossing points in (a), (b), (c) are taken to correspond the respective transition points for L — > oo (from 0). 

The universality class of the dynamic transition in Ising model under oscillating field has been studied extensively 
by investigating |5j the critical point and the cumulant value U* at the critical point, where the cumulant curves 
cross for different system sizes (L). In that case, of course, the variation of U (at any fixed L) is similar to that in 
the static Ising transitions (U = 2/3 well inside the ordered phase and U — ► well within the disordered phase). 
In fact, U* value in this oscillatory field case was found to be the same as that in the static case, indicating the 
same universality class 5J. We have studied the behavior of U in 2D. We observe different behavior in the field pulse 
induced magnetisation-reversal transition case. 

We observed Q two distinct behavior of the cumulant U. Typically, for low temperature and low pulse-duration 
region (see the inset in Fig. 2) of the magnetisation-reversal phase boundary, the cumulant crossing for different 
system sizes (L) occur at U* ~ 0.42 to 0.46 (see Fig. 8). We checked these results for both sequential and random 
updating. Specifically, for T = 0.5 and At = 5, (see Fig. 8c) for random updating we find the transition point 
value of hp ~ 2.6, to be smaller than the value (~ 1.9) for sequential updating. However, the value of U* at this 
transition point is very close to about 0.44. This indicates that updating rule does not affect the universality class 
(U* value), as long as the proper region of the phase boundary is considered. For relatively higher temperature and 
pulse-duration region of the phase boundary, the crossing of U for different L values occur for U* ~ + . This occurs 
both for sequential (Fig. 9a, b) and random (Fig. 9c) updating. It may be noted that the phase boundary changes 
with the updating rule, as the system relaxation time (which matches with the pulse width at the phase boundary) 
is different for sequential and random updating plj . 





FIG. 9. Behavior of U near the transition, driven by (a) h p at a fixed value of T (=2.0) and At (=5) with sequential updating, (b) T 
at a fixed value of h p (=0.5) and At (=10) with sequential updating, and (c) h p at a fixed value of T (=1.5) and At (=5) with random 
updating, for different L, averaged over 1000 to 6000 initial configurations. The insets show the typical behavior of the magnetisation m m 
at the time of withdrawal of the field pulse by varying (a) h p at a fixed T and At, for L = 50 and 400, (b) T at a fixed h p and At, for 
L = 50 and 200, (c) h p at a fixed T and At, for L = 50 and 200; m w = at the transition point (from Q). 

Note that in the low temperature and At regions, there seems to be significant finite size scaling of the transition 
(m w = 0) point (see the insets of Fig. 8a, b, c). In Fig. 8d, the finite-size scaling analysis of those data is presented. 
For the other cases, there seems to be no significant finite size effect on the transition point (cf. insets of Fig. 9a, b, 
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c) , indicative of a mean-field nature of the transition in this range. It may be noted that to compare the finite size 
effects, we normalise the parameters T or h p by their ranges required for full magnetisation reversal. In fact, this 
weak finite size effect for high T and At regions did not lead to any reasonable value for the fitting exponent in the 
scaling analysis. 

For the static transition of the pure two-dimensional Ising system, U* ~ 0.6107 0, For low temperature 
(and low At) regions of the magnetisation- reversal phase boundary, the observed values of U* (in the range 0.42 - 
0.46) are considerably lower than the above mentioned value for the static transition. There is not enough indication 
of finite-size effect in the U* value either (cf. Q). This suggests a new universality class in this range. Also, the 
finite-size scaling study for the effective transition points here (see Fig. 8d) gives a correlation length exponent value 
(y ~ 1.4) ?J larger than that of the static transition. This value of the exponent v compares well with that obtained 
by the data collapse (y — 2.0 ±0.5 in d = 2) ^4|. For comparatively higher temperatures (and high At), the U* ~ + 
at the crossing point. Such small value of the cumulant at the crossing point can hardly be imagined to be a finite-size 
effect; it seems unlikely that one would get here also the same universality class and U* value will eventually shoot up 
to U* ~ 0.44 (for larger system sizes), as for the other range of the transition. On the other hand, such low value of U* 
might indicate a very weak singularity, as indicated by the mean field calculations [Tflj mentioned in the introduction. 
In fact, even for the static transition, as the dimensionality increases, and the singularity becomes weaker (converging 
to mean field exponents) with increasing lattice dimension, the cumulant crossing point U* decreases (U* ~ 0.61 in 
d = 2 to U* ~ 0.44 in d = 4). We believe the mean field transition behavior here, as mentioned earlier, is even weaker 
in this dynamic case as reflected by the value U* ~ + , corresponding to a logarithmic singularity (as in eqns. ltTT|) 
and (fTTft S. 



V. SUMMARY AND CONCLUSIONS 



In this paper, we review all the aspects studied so far on the dynamic magnetization-reversal phenomena in the Ising 
model under a finite-duration external magnetic field competing with the existing order for T < T®. We have discussed 
the mean-field equation of motion and its solution, the nature of the MF phase boundary and the divergence of the 
susceptibility and relaxation time. The analytic solution of the linearized equation of motion also showed divergence 
of susceptibility and relaxation time on the MF phase boundary [T3| ■ The same transition has been studied by Monte 
Carlo simulations in both two and three dimensions. The obtained phase diagram is consistent with the classical 
nucleation theory. The nucleation process is initiated by the external magnetic field and depending on the strength 
of the field the system nucleates either through the growth of a single droplet or through the growth and subsequent 
coalescence of many droplets. For h p > h^sp the system belongs to the multi-droplet regime and the transition 
is continuous in nature; whereas for h p < hasp the system goes over to the single-droplet regime where transition 
is discontinuous. Expecting power law behaviour for both m w and £ in multi-droplet regime, the finite size scaling 
fits give the estimates of the critical exponents (3 and v for both d = 2 and 3. Unlike in the MF case, where the 
relaxation time t^ if shows a logarithmic divergence, tr in MC studies falls off exponentially away from m w — and 
the divergence in tr comes through the growth of the prefactor k with the system size. 

The universality class of the dynamic magnetisation-reversal transition in 2D, has been studied here using Monte 
Carlo simulations. Both sequential and random updating have been used. The phase boundary at any T (< T c °) is 
obtained first in the h p — At plane. They of course depend on the updating rule. The phase boundaries obtained 
compare well with the nucleation theory estimate h p \nAt — constant along the boundary 0,EJ. Extensive Monte- 
Carlo studies for the fluctuations in the order parameter \m w \ and internal energies etc. showed prominent divergences 
along the phase boundaries |6j. Fourth order cumulant (U) of the orderparameter distribution is also studied for 
different system sizes (upto L = 800) around the phase boundary region .7) . The crossing point of the cumulant (for 
different system sizes) gives the transition point and the value U* of the cumulant at the transition point indicates 
the universality class of the transition. In the low temperature and low pulse width range, the U* value is found 
to be around 0.44 (see Figs. 8a, b, c). The prominent discripancy with the U* value (~ 0.61) for the static 
transition in the same model in two dimensions indicates a new universality class for this dynamic transition. Indeed, 
the finite-size scaling analysis (Fig. 8d) suggests a different (larger) value of the correlation length exponent also. For 
comparatively higher temperatures and higher pulse widths, the U* values are very close to zero (see Fig. 9a, b, c), 
and the transitions here seem to fall in a mean-field-like weak-singularity universality class similar to that obtained 
earlier [l^, and indicated by eqns. (|llfl and (|17|) . Here, the finite size effects in the order parameter and the transition 
point are also observed to be comparatively weaker (see insets of Fig. 9). 
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